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Abstract: Bound states poles in scattering amplitudes are generated by the divergence of 
the perturbative series due to enhanced Coulomb scattering near thresholds. This suggests 
to organize bound state calculations according to an expansion in h, i.e., in the number 
of loops. I study QED and QCD bound states at lowest order in h, which are analogous 
to Born amplitudes. The absence of loops allows the use of retarded boundary conditions 
where particles only propagate forward in time, making a hamiltonian approach feasible. 
The instantaneous field is determined by the equations of motion separately for each 
Fock component of the bound state. The field equations are compatible with a linear 
potential as a homogeneous, non-perturbative solution. Stationarity of the action deter- 
mines the direction of the ensuing constant electric field in each Fock state. Applying this 
approach to relativistic qq and qqq states in QCD results in a bound state equation which 
provides a reasonable description of the spectrum, including linear Regge trajectories. The 
equal-time wave functions have unique Lorentz transformation properties, which ensure 
the correct dependence of the bound state energy on the center-of-mass momentum. The 
qqq potential is gauge covariant and confines the three quarks in a symmetric way. 



1. Introduction 



The perturbative expansion is a powerful tool in the analysis of scattering amplitudes for 
elementary fields (leptons, quarks and gauge bosons). Bound states (atoms and hadrons) 
appear as poles in scattering amplitudes which arise from the divergence of the perturbation 
series. The study of bound states and their interactions thus requires summing an infinite 
set of Feynman diagrams. This may be done in the Bethe-Salpeter framework [1] . However, 
the choice of diagrams to be included in a first approximation, as well as the systematic 
ordering of the inclusion of the remaining contributions, allows a wide range of approaches, 
particularly for relativistic bound states. 

Here I shall argue that an h (or equivalently, loop) expansion may serve as a guiding 
principle in bound state calculations. In the ?i — t- limit the phase exp{iS/h) in the path 
integral oscillates rapidly, causing the (bosonic) fields to approach their classical values, 
corresponding to a stationary action S. In the absence of field fluctuations there are no 
loop integrals - thus the Feynman diagrams of lowest order in h are Born terms. Since 
each loop brings a factor hg'^ the perturbative and loop expansions are equivalent for a 
given Green function'^. Bound state wave functions on the other hand contain all powers 
of g, yet may have no loop contributions, thus decoupling the h and g expansions. This 
is illustrated by the way in which the Schrodinger and Dirac equations emerge from QED 
perturbation theory. 

The Schrodinger equation for non-relativistic (NR) atoms like muonium {e~ is 
generated by the sum of ladder diagrams shown in Fig. 1. The factor a = e'^/An suppression 



Figure 1: Sum of ladder diagrams in e /i+ scattering which generate a non-relativistic bound 
state pole. 

of each photon exchange is compensated when the propagators in the corresponding loop 
integral are nearly on-shell, which brings an inverse factor of the relative velocity v between 
the fermions. Thus for f ~ a all diagrams are effectively of the same order, causing the 
series to diverge regardless of how small a is. Due to the NR motion Coulomb photon 
exchange dominates and transfers only 3-momentum. Energy transfer is suppressed by 
a, quenching the loop integrals in Fig. 1 at leading order. Hence the dynamics reduces 
to scattering from a fixed a/r Coulomb potential and is described by the Schrodinger 
equation. The repeated scattering from the potential brings powers of a without loops. 

The Dirac equation for a relativistic electron may be derived in a similar way [3], by 
including all crossed photon exchanges as indicated in Fig. 2. In the limit of a large muon 
mass the full sum of crossed and uncrossed photon exchange diagrams reduces at each 
order in a to the scattering of the electron in a fixed a/r Coulomb potential described by 



^The equivalence between the number of loops and powers of h is actually non-trivial, see [2] 
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Figure 2: Sum of all uncrossed and crossed photon exchange diagrams which generate Dirac bound 
states in the limit where one of the fermion masses is large. 

the Dirac equation^. Hence the leading contribution of the sum of multiloop diagrams in 
Fig. 2 to the bound state poles arises from classical Coulomb scattering and in this sense 
is of lowest order in h. 

The symmetries of the lagrangian are preserved at each order in h. Thus boost invari- 
ance requires that the bound state energy E must be related to its CM momentum k as 
E = \/m2" + fc^. Wave functions defined at equal time in different frames must also be 
related, but that relation is dynamical since the concept of equal time is frame dependent. 
I find an example of this in the present approach. 

Scattering amplitudes at the Born level are insensitive to the ie prescription of the 
propagators. The relative magnitudes of the time-ordered contributions to these amplitudes 
on the other hand do depend on the prescription. Bound states evaluated at lowest order 
in h are analogous to Born terms. I show that the energies of Dirac bound states indeed 
are independent of the ie prescription, whereas the equal-time wave functions do depend 
on the prescription. The standard single-electron Dirac wave function corresponds to using 
a retarded electron propagator, which only allows forward propagation in time. The use 
of a Feynman propagator would give Fock states having arbitrarily many electron-positron 
pairs, arising from backward propagation of the electron in time (Z- graphs). The use of 
retarded propagators thus essentially simplifies bound state analyses at lowest order in h. 

Gauge theory lagrangians have no time derivative of ^4*^ (unless this is introduced via 
the gauge fixing). In the ?i — ?• limit is determined at each instant of time by the 
positions of the charged constituents. This implies that Fock states with constituents at 
different positions have distinct fields. The field equations allow a linear A^ potential as 
a homogeneous (non-perturbative) solution. For the action to be stationary the orientation 
of the potential needs to be correlated with the locations of the charges in each Fock 
component. Consequently a linear potential is allowed in a self-consistent semi-classical 
treatment, maintaining Poincare invariance. 

In the next section I outline the main results, which are then discussed in more detail 
in the following sections. 

2. Outline 

2.1 QED atoms and the linear potential option 

I use a hamiltonian approach to field theory bound states in the semi- classical limit, i.e., 
at lowest order in h (no loops). I first illustrate the procedure for ordinary QED atoms 

■^In a Bethe-Salpeter approach the crossed photon exchanges require kernels of unlimited orders in a. 
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(muonium, e^/i^) in their rest frame, whose binding is known to be dominated by the 
Coulomb potential. This example reminds us that the instantaneous potential 
determined by the equations of motion is different for each Fock state component. It also 
clarifies the distinction between the A^ field for a given Fock state and the A^ field that 
would be measured by an external probe, the latter being given by a superposition of all 
contributing Fock states. 

The determination of ^4*^ from the QED field equation for the Fock component where 
the electron is located at xi and the muon at X2, 

-Vl A^{x; xi,X2) = e6^{x - xi) - e5^{x - X2) (2.1) 

is unique only up to a homogenous solution. In particular, a contribution linear in x may 
be added, 

A\x- x^,X2)=K^i-x + -^ i\Zr^\ - uT^l = + (2.2) 

47r \|a; — a;i| \x — X2\ J 

where the magnitude A and the direction t of the linear term are free parameters. The 
gauge field (2.2) gives an instantaneous contribution to the action. 



1 

47r \xi — X2\ 
(2.3) 



Here I suppressed the infinite self-energy contribution of the Coulomb potential, which is 
independent of xi and X2 and thus irrelevant for bound state evolution. The first term in 
(2.3) is likewise infinite and irrelevant provided A is a universal constant, A ^ A{xi,X2). 
The interference term between the linear and Coulomb contributions is finite and may be 
evaluated by taking the integration region to be a sphere whose radius R tends to infinity. 
Using ^'^A^-^ = and Stokes' theorem, 

[d'xV- [{VAUAi^i] = lim / dS ■ (V^L)^Li = \eA^ ^ (^1 " ^2) (2.4) 

where terms of O (l/i?^) in the Coulomb potential could be neglected, 

1 1 r X ■ Xi~\ 

= — IH , - 

\x-Xi\^^^^^ Rl i?2 J|3.|=^ 

Stationarity of the action (2.3) wrt. variations in the direction of the linear potential 
requires £ to be parallel to Xi — X2- Choosing the sign to give an attractive (rather than 
repulsive) potential leads to 

,n/ . . Xt — X-? e / 1 1 \ , , 

A^{x;xi,X2) = ^-j-x + —[- -] 2.6 

\Xi — X2\ AtT WX — Xi\ \X — X2\) 

with A a universal constant. The linear potential allows to incorporate non-perturbative 
features of the theory at the semi-classical level. The fermion interaction with the field 
gives 

e \A^[xx\xx,X2) - J^{x2\xx,X2)\ = eA^ \xi - X2\ - 2—- r (2.7) 

47r \Xi — £C2 



+ 0(4) (2-5) 
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where I again omitted the infinite, xi- and a;2-independent Coulomb self-interactions. 
The Coulomb potential is counted twice since both the electron and the muon inter- 
act. However, according to (2.3) the phase originating from the gauge part of the action 
exp(— I f dt f d^xF^yF'^'^), which describes the energy stored in the Coulomb field, reduces 
the Coulomb potential to its physical value, Vc = —a/r. 

The experimental fact that QED does not confine electric charge compels us to choose 
^QED = and thus recover the standard description of muonium bound only by the 
Coulomb potential. In Sections 5 and 6 I find corresponding linear potential solutions of 
the QCD equations of motion for baryons and mesons. Choosing Aqcd 7^ amounts to a 
novel boundary condition for perturbation theory, motivated by the confinement of color 
in QCD. 

A linear potential extending to infinity would be unacceptable on physical grounds. 
However, A^{x] xi,X2) of (2.6) is the field of only a single Fock state. The field that 
would be measured by an external probe is a superposition of the contributions from all 
Fock components. As seen from (2.6), Fock components with opposite separations xi — X2 
give opposite contributions to A9 which cancel if the corresponding Fock probabilities are 
the same. Hence the electric field vanishes outside the bound state, where all Fock states 
contribute coherently. 

The linear potential contributes at O (e) since A is non-perturbative, see (2.7). Hence 
in a perturbative sense it dominates the single (Coulomb and transverse) photon exchanges 
which are of O (e^) . This allows the analog of a "Born term" for gauge theory bound states, 
which is exact at lowest order in h (no loops) and at first order in the coupling. Lorentz 
invariance must then be preserved. In particular, the bound state energy E should have 
the correct dependence, E = + k^, on the CM momentum k. The bound state 

equation that we derive here for wave functions at equal time in all frames in fact has this 
non-trivial and unique property. By contrast, in QED atoms the Coulomb potential 
dominates transverse photon exchange only in the rest frame. The boosted muonium atom 
gets a leading O (e^) contribution from \efi'y) Fock states with a transverse photon which 
must be included to get the correct k dependence [4]. 

2.2 Dirac equation in QED — retarded boundary conditions 

Particle production can no longer be neglected when the dynamics is relativistic. An |e/i) 
state at t = will thus with time develop into states with additional electron and muon 
pairs if the interaction potential is commensurate with the particle masses. Relativistic 
bound states (defined at equal time of the constituents) therefore have Fock states with 
arbitrarily many particles. 

The Dirac wave function ip{x) of an electron bound in an external potential is relativis- 
tic yet describes the spatial distribution of a single (positive or negative energy) electron. 
Understanding the apparent absence of multi-particle Fock states in the Dirac wave func- 
tion turns out to be useful for formulating a relativistic hamiltonian description of e~ 
atoms and hadrons. 

The Green function G{p^,p) of an electron in a static (time-independent) A^ potential 
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(Fig. 3) satisfies 

G{p^,p) = S + SKG (2.8) 
where S is the electron propagator and K the kernel for a single interaction. The 

p\0 p\p S K K S 
S > <! > <■ > » ^ > ^ 

G(p^,p)= + Cp + k2 <^j^^ + S S ^ +... 



Figure 3: Electron scattering from a static external potential. The energy of the electron does 
not change during the scattering. The initial and final electron momenta are denoted {p^,0) and 
P = {P°,P)- 

component of the electron's 4-momentum is conserved during the scattering since the static 
source only transfers 3-momentum. I take the initial and final 3-momenta of the electron 
to be and p, respectively. 

The Green function has poles at the bound state energies Epi, 



G(Ap) = 9^ + ... (2.9) 



R{Er,P 
P^-Er 

The pole residue R{Eji,p) satisfies R = SKR, or explicitly (with p^ = Er) 

R{Er,p) = — / (-ie)j^A\k)R{ER,p-k) (2.10) 

Multiplying hy p — m and Fourier transforming to {p^, £c)-space we find that the residue R 
satisfies the Dirac equation, 

iV • 7 + e7°^°(ic) + ml R{Er, x) = Er^j^ R{Er, x) (2.11) 

In order to display the equal-time Fock states of the bound state given by the Dirac 
equation we need to time-order the interactions. At O (e^) the first (second) diagram on the 
rhs. of Fig. 4 corresponds to the intermediate electron having positive (negative) energy, 



Ei = zizy kf + m'^. According to the Feynman ie prescription this electron propagates 
forward (backward) in time, corresponding to an intermediate \e~) {\e~e~^e~)) Fock state. 
At higher orders in e further time orderings contribute. Consequently the bound state has 
Fock components with arbitrarily many e~^e~ pairs. The creation and destruction of the 
various Fock components balance to create a bound state which is stationary in time. 

The Dirac wave function il^{x) describes only a single electron, not the details of the 
multiparticle bound state dynamics. As seen above, the bound state energies = Er 
given by the poles of the Green function (2.9) nevertheless are eigenvalues of the (time 
independent) Dirac equation (2.11) satisfied by the pole residues R{Er,x). This single 
particle Dirac wave function may be given the following hamiltonian interpretation. 

The static potential conserves the energy component > of the electron momentum. 
Hence the covariant Feynman diagrams in Fig. 3 which build the Green function G{p^,p) 
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E.>0 E.<0 




Figure 4: Electron scattering from a static A" potential at O (e^). The covariant (p*^-conserving) 
diagram on the Ihs. splits under time-ordering into the two diagrams on the rhs. The interme- 
diate electron has 3-momentum ki and energy Ei = -ii\J + m^. The negative energy electron 
propagates backward in time, implying the creation (annihilation) of an e^e~ pair at ti (^2)- 



do not depend on the Feynman ie prescription at the negative energy pole of the electron 



propagator, p*^ = — y fc? -|- m1 + ie. In particular, the bound state energies Er of the Green 
function (2.9) will not change if instead of the Feynman propagator of the electron we use 
the retarded one. 



Sr{p^,p) = i—^ — ^—rr4 — r (2.12) 



where Ep = yjp^ + rrig. However, this changes the time development of the scattering so 
that the electron always propagates forward in time, 

SR{t,p) = ^ [{Ep^^ -P-1 + me)e-^^-* + (i?p7° + P • 7 " meje^^''*] (2.13) 

Consequently the second (Z-)diagram on the rhs. of Fig. 4 is absent, while the first dia- 
gram refers to the (positive or negative energy) electron moving forward in time. There is 
only a single electron Fock state and the Dirac wave function describes the distribution of 
this electron. This single particle picture of the time development obtained using retarded 
boundary conditions does not correspond to the physical Fock state structure. Further- 
more, the argument that the bound state energies {Epi in (2.9)) are independent of the 
boundary condition only holds in the absence of loops, i.e., in the ^ — ?• limit^. 

In the absence of backward propagation the electron becomes localized in space as the 
propagation time vanishes, 

lim SRit,x) = -i^S^ix) (2.14) 

which makes a hamiltonian approach feasible. In the operator formalism the retarded 
boundary condition can be implemented using the "retarded vacuum" , 

|0)i? = iV-^n<AlO) (2.15) 



^Loops with retarded propagators vanish since all poles are on the same side of the energy axis. 
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where the product is over all momenta p and helicities A. The (infinite) normahzation 
factor is fixed by i?,(0|0)R = 1. The retarded vacuum satisfies 

bp,x\0)R = dljO)R = and hence V(x)|0)ij = (2.16) 

where 'ip{x) is the free (interaction picture) fermion field. Consequently the retarded prop- 
agator (2.13) is given by the standard operator matrix element in the retarded vacuum, 

Sr{x -y) = r{Q\ T[^{x)i^{y)\ \^)r (2.17) 

The negative energy contribution to the propagator arises from the d)d term, which repre- 
sents the removal of a positive energy antifermion from \^)r. The interaction hamiltonian 
annihilates the retarded vacuum, 

Hi{t)\Q)R = e j d^xA^{x)i)\t,x)'>lj{t,x)\{))R = {) (2.18) 

which ensures the absence of particle production. In Section 4 I derive the Dirac equation 
using this hamiltonian approach. 

2.3 Relativistic QCD bound states 

The instantaneity of the potential is a consequence of the absence of a d^A^ term in 
the lagrangian. The dominance of A^ over the transverse A gauge field components in a 
perturbative expansion holds for non-relativistic motion (atoms in their rest frame) and 
more generally with a non-perturbative linear ^ potential as in (2.7). At lowest order 
in h the use of the retarded vacuum (2.15) as boundary condition maintains bound state 
energies while allowing a valence wave function, as discussed above for the Dirac equation. 

In Section 5 I find a QCD meson solution using a hamiltonian approach in the retarded 
vacuum. In a specific gauge the linear potential (2.6) is present only in the commuting 
gauge fields {A^ and A^) of SU3. The equation is expected to be exact at leading orders in 
h and the coupling 17, giving the equivalent of a Born term for bound states. The qq color 
singlet wave function X}f{x), where k is the CM momentum of the bound state, a and /3 
are Dirac indices and x is the spatial separation of the quark pair at an instant of time, 
satisfies the bound state equation 

-iV • [a, Xk{x)\ + \k ■ {a, Xk{x)] + nii-f^Xkix) - Xkixh°m2 = {E - V)xk{x) (2.19) 

Here mi (rfi2) is the mass of the q (q) constituent, o; = 7*^7 is a standard Dirac matrix and 
V{x) = gA'^\x\ is the linear potential. This equation is similar in form to the one proposed 
by Breit [5] already in 1929. Here the potential is purely linear since the derivation is valid 
only at leading order in the gauge coupling g. The properties of this equation (with fc = 
and mi = 711,2) was studied phenomenologically [6] for a linear + Coulomb potential (see 
also [7] and references therein) . It was also derived using a stationary phase approximation 
with retarded boundary conditions [8] . 

Due to the underlying gauge invariance it is perhaps not surprising that the energy 
E of the bound state appears in the canonical form E — eAP in (2.19), while the CM 
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momentum k is kinematical since A = 0. The radial wave functions have singularities at 
r = and at -E — V{r) = 0. Requiring local normalizability at these points gave quantized 
energy levels and a phenomenologically relevant spectrum, including asymptotically linear 
Regge trajectories [6]. The bound state equation reduces to the Schrodinger equation in 
the non-relativistic limit {V <^ E)^ in which case the singular point E — V{r) = moves 
to r = CO. 

At large quark separations r the linear potential V{r) dominates the rhs. of (2.19), 
and can be balanced only by a large derivative term on the Ihs. Hence the qq wave function 
is rapidly oscillating at large r, x ~ exp(zgA^r^/4), and has an r-independent probability 
density. Since this wave function was derived using retarded boundary conditions it in- 
directly (via its negative energy components) describes the virtual "sea" quarks obtained 
with physical boundary conditions. The normalization oc r of the wave function (2.19) 
plausibly reflects the quark pairs created by the large potential energy V{r) oc r. This 
question deserves further study. 

Since the wave function Xk{x) describes constituents at equal time for all CM momenta 
k it is not explicitly Lorentz covariant. Nevertheless, if the bound state equation (2.19) is 
accurate to lowest order in h and g the bound state energy must have the correct dependence 
on the CM momentum, E{k) = \/ + M"^. Remarkably, this turns out to be the case [9]. 
The fc-dependence of the wave function is explicit in 1+1 dimensions. In 3+1 dimensions 
Xfc(iE) with a; II can be expressed in terms of the fc = wave function. This also holds 
for the first derivative of Xfc n^^t- x k, which (numerically) allows to determine Xk{x) 
for all X from the = solution. 

The fc-dependence of the wave function Xk ^-nd of its energy eigenvalues E{k) are 
found as follows [9]. For fc = (0,0, k) along the z-axis, the bound state equation (2.19) for 
Xa:(0,0,z) and its first derivatives 3jXfc(0, 0, z), (i = x,y) become independent of k when 
the coordinate z is expressed in terms of the "boost invariant" variable s defined by 

In the non-relativistic regime, V E, this is just the Lorentz contraction expected for 
distances measured at equal time in every frame. The fact that the transformation is 
governed by the canonical energy E — V{z) means that the contraction rate for a relativistic 
wave function is z-dependent. The Lorentz covariance of (2.19) requires furthermore that 
the square of the canonical four-momentum P = (£^ — V, 0, 0, k) is frame independent when 
expressed as a function of s. For a linear potential V{z) = ^A^l^l the relation (2.20) 
integrates (taking z > 0) to 

s = \z{E-\V) (2.21) 
With E"^ — k'^ = M"^ we then have indeed that 

= {E- V f -k^ = - Agh?s (2.22) 

is independent of k. This seeming coincidence only holds for a linear potential. 

For baryons I find a corresponding qqq bound state equation and a gauge where the 
wave function is fully antisymmetric in color. Stationarity of the action implies that the 
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color electric fields VA3 and V^g point in separate directions, thus confining all three 
quarks. The relative strength ylg/Ag depends on the positions xi,X2,X3 of the quarks in 
such a way that the potential is the same for all permutations of quark colors. This allows 
(in a specific gauge) the simple color structure eABC of the wave function, as assumed in 
the quark model. 



3. Non-Relativistic QED Bound States 

Some aspects of the methods sketched above may be illustrated in the familiar case of 
non-relativistic (NR) QED atoms, specifically muonium (e~//^). I use an equal-time, 
hamiltonian field theory approach which is close in spirit to the usual treatment in quantum 
mechanics. However, rather than reducing the problem to relative motion in a fixed central 
potential I determine the instantaneous field by taking matrix elements of the QED 
operator equations of motion, 

d^F^^'^ix) - e ^ V5i(x)7>i(x) = (EOM) (3.1) 

i=e,fi 

The resulting depends on the matrix element, and specifically on the positions of the 
electron (a^i) and muon {x2)- Hence each Fock component of the bound state is assigned a 
distinct potential A^(x; xi, X2). This is made possible by the instantaneity of A^, which in 
turn is a consequence of the absence of a time derivative BqA^ in gauge theory lagrangians. 
The Coulomb gauge condition (V • A = 0) maintains the instantaneity of ^ and in this 
sense is preferable to covariant gauge fixing {dfj,A'^ = 0). 

In the NR limit we may ignore fermion pair production (I address the relativistic effects 
in Section 4). The dominance of the A^ potential^ and its instantaneous nature then implies 
that only two-body \e~ n'^) Fock states contribute at lowest O (e^) in the coupling. 

I parametrize the e~fi'^ bound state at time t = as 

\E,t = o) = j d'y^Sy^ i^l{t = 0, yi)x{yi,y2)Mt = 0, 2/2)10) (3.2) 

where the wave function x{yi->y2) is a 4 x 4 matrix in Dirac space^. The free fermion 
operators of the Interaction Picture usual 



V'(x) = / (^^^^ E ^P' ^)e"''"^,A + v{p, A)e^^-4^,] (3.3) 



and satisfy the canonical anticommutation relation 

{Va(t, x), i,\{t, x')] = 5\x - x')5^p (3.4) 

''This holds only in the atomic rest frame. For atoms in relativistic CM motion also I//7) Fock states 
with a transverse photon contribute at lowest order [4] . 

use a relativistic notation here and later reduce the wave function to its leading NR components. 
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The matrix elements of the EOM (3.1) for each |e fi'^) Fock component of the bound 
state (3.2) at t = should vanish, 

(0|V^]^^(0, X2)^Pea{0, xi) (EOM) \E , 0) = (3.5) 

Since the Fock states do not contain physical (transverse) photons only the classical (in- 
stantaneous) field survives in F^'^{x) at lowest order in the coupling e. For ;/ = in 
(3.1) the constraint (3.5) is 

Xap{xi,X2) [-V^A\x)\ = e ^(0|v],^(0,a;2)V'ea(0,a;i)vl(0,a;)V^,(0,a;)|S,0) 

i=e,ii 

= e [8^{x - xi) - 5'^(x - X2)\ Xap{xi,X2) (3.6) 
where I used (neglecting multi-particle states in the NR limit, hence ipe{x)\Q) — ?• 0) 

V'e(0, x)^M^, x)\E, 0) = Vl(0, x) j d^y^xix, 2/2)V'm(0, ^2)10) (3.7) 

The standard solution is 

A\x; x,,X2) = ^ f - r^—) (3.8) 

47r\|ic — a^il \x — X2\J 

The interpretation of this result differs from the case where A^{x) is regarded as a fixed 
external Coulomb potential centered at x, which is sampled by the charged constituents 
according to their positions xi and X2- Now there is no external field but rather a bound 
state (3.2) which is a superposition of Fock states. The gauge field A^ is constrained 
for each Fock component and each instant of time by the QED equation of motion. At 
large \x\ the Coulomb field (3.8) vanishes faster than l/\x\ as appropriate for a neutral 
state. Furthermore, (3.8) is the field of a single Fock component. A measurement of A^ 
by an external probe would involve a sum over all Fock contributions weighted by their 
probabilities. This ensures, e.g., the cancellation of the dipole field outside 5- wave bound 
states. 

For v = j (= 1, 2, 3) in the EOM (3.1) the A^ field does not contribute in the matrix 
element (3.5), 

Xap{xi,X2) ■ = e5^{x - iCi)(7°7^x)a/3 - e5^{x - X2){xi^l^)ap = (3.9) 

This is satisfied at lowest order in the coupling e for NR QED atoms since the space 
components j of the constituent currents are suppressed compared to their charge densities 
by a factor p/m ~ a, where p is the relative momentum and m the mass scale of the 
constituents. In analogy with the scaling of the upper and lower components of a Dirac 
wave function I expect (and verify below) that the wave function x{^ii^2) has the Dirac 
structure, in 2 x 2 block form. 




X= ^'M , <; Xii, X22 = 0{p/m) (3.10) 

\ X21 X22 / 



0{l) 

1 {p/m 
O {p^/m'^) 
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For the dominant Fock components described by xi2 the rhs. of the EOM (3.9) is thus 
suppressed by p/m due to the off-diagonal nature of 7*^7-'. The EOM (3.6) for the other 
components Xij is of higher order in due to the suppression of the corresponding wave 
functions. 

Having determined the field (3.8) for each Fock component and verified that the 
EOM is consistent with A = we may impose stationarity in time on the Fock amplitudes, 



4>al3{t; Xi,X2) = {0\i)j^p{t,X2)i)ea{t,Xi)\E,t) = 6 *^Va/3(i = 0;Xi,X2) 



(3.11) 



where (papit = 0;a;i,a;2) = Xai3{xi,X2) follows from the anticommutation relation (3.4). 

The time dependence of the fermion fields is explicit in (3.3) while that of the state 
\E,t) is given by the interaction Hamiltonian 



Hi{t) = e J (fxA'^ix) i)]{t,x)i)f{t,x 



(3.12) 



/=e.M 



Hj actually appears with a factor ^ in Coulomb gauge quantization (see, e.g., Appendix 
A of [10]). Equivalently, and perhaps more intuitively, we may consider (3.12) as giving 
the interaction energy of the electron and muon with the Coulomb field, but then need to 
consider also the energy stored in the field, since it depends on the Fock state. 



exp 



idt / d'x{-\F^,F^'') 



exp 



exp 



idt 
-idt 



j d^xA^{x)V'^A^{x) 

exp{-idtEA) (3.13) 



1 



AtT \xi — X2\ 



where I discarded an a;i, a;2-independent infinite phase (self-energy). This contribution 
cancels half of the e,/x interaction energies arising from (3.12). 

The stationarity condition for the bound state wave function (3.11) is then (at t = 0), 



y/3{0;xi,X2) 
dt 



(0|i- 



Vea(0,a;i)|i?,0) + i^lp{^,X2] 



d'ljjeaiO,Xi] 



dt 



+ {0\^Plg{0,X2)^Pea{0,Xl)[HJ{0) + Ea]\E,0) = E (l)a^{0;xi,X2) 



\E,0) 
(3.14) 



Using (3.7) and (3.4) gives 

{0\^lJl^{0,X2)^PeaiO,Xl)HIiO)\E,0) = [eA'^ixi) - eA'>{x2)] xixi,X2) 



1 



An \ xi — X2 



-X{Xl,X2) 



(3.15) 



where I again discarded the infinite self-energy terms. Due to the electron interacting with 
the field of the muon and vice versa the potential energy is twice the physical value. In the 
bound state equation (3.14) this is reduced by the field energy Ea of (3.13) to the standard 
Coulomb potential 



V{xi,X2) 



1 

AtT \xi — X2\ 



(3.16) 
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Altogether the stationarity condition (3.14) gives, 



7°(-zVi-7 + me)x(a;i,a32)-x(a^i,a^2)7°(^V2-7 + m^) = [E -V{xi, X2)]x{xi, X2) (3.17) 

where V2 = 8/8x2 operates to the left. This is identical in form with the relativistic bound 
state equation that I derive in Section 5, but presently valid only in the non-relativistic 
limit where particle production and the transverse components of the gauge field A may 
be ignored. 

The reduction of (3.17) to the non-relativistic Schrodinger equation is analogous to 
the standard reduction of the Dirac equation. In the CM, xixi,X2) — )■ xi^i ~ ^2) = x{^)- 
Ignoring the O (p^/m^) component X21 in (3.10) the bound state equation (3.17) is in 
component form 



iVxi2 • o- -ia- ■ Vx22 + « Vxii • <^ | _ / 2m^ xii (-^b - ^)xi2 

-icr • Vxii + « Vx22 • o" -io-'Vxu / V 2meX22 



(3.18) 



where on the rhs. I defined the O {p^/m^ binding energy Eb through E = rrif. + + E^, 
and approximated E — V — me + ~ 2m^, E — V + rrie — ~ 2me. Substituting 

Xii = 7, Vxi2 o- X22 = 7, — V • (^Xi2 (3.19) 



in (3.18) gives the Schrodinger equation for xi2, 

+ v) X12 = EbXi2 (3.20) 



2me 2m^ 



As expected the Schrodinger equation is independent of the 2x2 matrix structure of xi2; 
and hence of the electron and muon spin orientations. 



4. The Dirac Equation 

The electron equation of motion, 

{i0 - e4 - me)'ip{x) = (4.1) 

follows directly from the QED lagrangian and has precisely the form of the Dirac equation. 
However, (4.1) is an exact operator relation whereas the Dirac equation gives an approx- 
imate c- numbered amplitude ■iIj{x) for electron scattering in an external gauge field. The 
Dirac equation may be derived from the field theoretic Bethe-Salpeter equation by consid- 
ering e~ scattering in the limit where the muon mass tends to infinity [3]. Somewhat 
surprisingly it turns out that the interaction kernel cannot be restricted to any finite or- 
der in the coupling e but must include an unlimited number of crossed photon exchanges, 
as also seen from Fig. 2. On the other hand, the Dirac equation does not include loop 
corrections on the electron propagator or vertices. 

As already discussed in Section 2.2, and pictorially shown in Fig. 4, the interaction 
hamiltonian (3.12) creates fermion pairs in the perturbative vacuum. Nevertheless, the 
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Dirac (c-numbered) wave function 'tp{t, x) at a given time depends on the coordinate a; of a 
single (positive or negative energy) electron. This may be understood by considering elec- 
tron scattering in the retarded vacuum (2.15), \0)r = iV~^ Yip \ '^p aI*-*)' where the product 
extends over all positron momenta p and helicities A. With retarded boundary conditions 
both positive and negative energy electrons propagate forward in time, Z-diagrams are 
absent and according to (2.18) there is no particle production, 

Hi\0)r = (4.2) 

In a covariant perturbative description of electron scattering from a static potential 
A^{x) the electron energy is unchanged by the scattering. Hence, for > and in the 
absence of loop integrals, the ie prescription at the negative energy pole p^ = —\/p'^ + ml 
of the electron propagator S{p^,p) is irrelevant. Consequently the bound state energies, 
i.e., the locations of the poles in the S-matrix, are the same for Feynman |0) and retarded 
\0)r boundary conditions. On the other hand, the Fourier transform p*^ — t- i is sensitive to 
the ie prescription at the negative energy pole. The Dirac wave function ip[t,x) describes 
the single electron bound state obtained with retarded boundary conditions, thus hiding 
its true multi-pair Fock state structure. 

It is straightforward to derive the standard Dirac equation for an electron bound by an 
external potential AP{x) using the hamiltonian method of Section 3 in the \^)r vacuum. 
The bound state at f = is parametrized in terms of its 4-component wave function '^{x) 
as 



\E,t = d)= / (fx'ilj\t = Q,x)ip{x)\{))R 

Sp 

A 



^2^^3 2Ep ^ K^^' ^)'^(P)^i,AlO)^ + ^^(-P' A)^(p)d_p,A|0)K] (4.3) 



The negative energy components of (p{p) describe a state where d_p x has removed a 
positive energy antifermion from \0)r. For the bound state to be stationary in time each 
Fock state amplitude (j){t, x) must satisfy 

ct){t, x) = R{0\^P{t, x)\E, t) = e-^^V(0, x) (4.4) 

where 0(0, a;) = ip{x) follows from ^ipait, x) , 'ip'p{t, x')^ = 6^{x — x') 6ai3 and iP{x)\0)r = 0. 

The time dependence of \E,t) is given by the interaction Hamiltonian (3.12). The 
stationarity requirement for the bound state at t = is then 

.d^{^ = R{0\^^^f^\E, 0) + R{0\m x)Hi\E, 0) = i?c/.(0, x) (4.5) 

The interaction picture fields satisfy 

i^^^^^ =/(-iV-7 + me)V'(t,a;) (4.6) 

and, making use of i{(0|V'^(t, x) = 0, 

r(0|V'(0, x)Hi\E, 0) = eA^xMx) (4.7) 



-13- 



Using these relations in (4.5) gives the Dirac equation for the wave function ip{x) of a 
bound state of energy E in the external potential A^(x), 

{-iV ■ 7 + e7°A°(a;) + me)(p{x) = E-f^ip{x) (4.8) 

The time development of fermion states may equivalently be derived using functional 
integral methods [11]. The retarded vacuum \0)r of (2.15) is then seen to correspond to a 
vacuum without a Dirac sea. I briefly recall this in the Appendix. 

5. Meson bound state equation in QCD 

I now apply the hamiltonian method demonstrated for non-relativistic QED atoms in 
Section 3 to relativistic color singlet ud mesons, using the standard QCD lagrangian 

Cqcd = -iFrF^, + ^ ijf{i0 - gi,T%s - m/)^f 

f=u,d 

Fr = d^A'^ - d^A^, - gfabcA^.A'^c (5-1) 

The explicit description of particle production in the perturbative vacuum is avoided by 
using the retarded vacuum, shown in Section 4 to give the Dirac equation and correct 
bound state energies at lowest order in h. The filling of all u and d states now includes a 
product also over colors A, 

\0)r = N-' n d^\p,X)bf{p,Xm (5.2) 

ensuring that the interaction picture fields satisfy 

i'^it,y)\0)R = tl:f{t,y)\0)R = (5.3) 
The filling of all u states may be equivalently expressed in coordinate space, 

H d\{p,X)\0)^lli^A{t,x)\0) (5.4) 

p,A,A x,A 

A local SU(3) gauge transformation U transforms ipAi^) — ^ U ab{x)''P b{x) and hence 

V'l "02 V'S V'l i'2 1p3 ^ £ABC U1AU2BU3C = tpl V'2 1p3 det U = Vl V'2 -03 (5.5) 
A,B,C 

where the antisymmetric tensor sabc arose from anticommuting the fermion operators. 
Hence the retarded vacuum (5.2) is a color singlet. It is also Lorentz invariant, since 
U\d'^{p, A) C/y^^ = d'^{A.p, A) only amounts to a reordering of the momenta. 
The ud bound state at t = is similarly to (3.2) expressed as 

t = 0) = / d^'y^d'^y^ V«^(t = 0, yi)x^''{yi,y2)i^d{t = 0, ^2)10)^ (5.6) 
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The wave function describing the state (5.6) is gauge dependent. For the state to be 
invariant under time-independent gauge transformations 'tp{t,x) — t- U{x)ip{t,x) we need 



xiyi,y2) ^ Uiy,)xiyi,y2)u\y2) (5.7) 

The stationarity condition for the bound state wave function corresponding to (3.11) is 

ct>'^^{t;xi,X2) ^ Rm!j'^{t,X2)i;^^{t,xi)\E,t) = e-'^'cl^'^^it = 0-,xi,x^ (5.8) 

where (pa^it = 0; cci, X2) = Xa/si^i, X2) follows from (5.3) and the anticommutation relation 

[^t{t,x),^f{t,x')] = 5\x - x')5^p5^'' (5.9) 

As in the non-relativistic case (3.5), matrix elements of the QCD operator equations of 
motion should be satisfied to leading order in g for each Fock state, 

R{0|V'f;(t,a;2)V'S.(i,=«i)[9Mi^a^" + 5/a6cFr^^-5 E V^/7X^''V'f]|i?,t) = (5.10) 

f=u,d 

From (5.3) and (5.6) we then get the following constraints on the gauge fields at t = 0: 

X^^(a;i,a:2) [d^Fr + gfabcF^'Al] = g ^^(^ _ Vt ^2) 

- g6\x - X2)x''^{xi,X2h''^''T^'' (5.11) 

The sources generate fields of O ((7), whose interactions contribute to bound state evolution 
at order {e.g., single gluon exchange between the quarks). However, the EOM (5.11) 
allows also an O (5^) instantaneous field 

Al{x)=Alia-x (5.12) 

as a homogenous solution (V^yl[J = 0) which contributes to bound state evolution at O (g). 
I shall consider only terms of this order, and thus neglect the O (5^) contributions from 
the vector potentials Aa {j = 1, 2, 3). As an ansatz for a bound state solution I assume the 
existence of a gauge where the wave function is diagonal in color, 

x^''iyi,y2) = s^''xiyi,y2) (5.i3) 

In order that the quark interactions with the gauge fields do not generate color non-diagonal 
Fock states only the commuting a = 3, 8 components of the A'^ gauge fields (5.12) can be 
non- vanishing in the gauge (5.13). 

The action gets an O (g) contribution from the interference between the O {g^^ fields 
(5.12) and the O (g) fields ^0=33 generated by the sources in the EOM (5.11). Due to 
/a38 = the gluon self-interaction term does not contribute at lowest order. With z/ = 
and quarks of color C the EOM constraint becomes 

-V^A'iix) = gT^^' [6\x - x,) - 6\x - X2)] (5.14) 
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(no sum over the quark color C). With the homogeneous solution (5.12) we have then for 
a = 3, 8 the instantaneous potential 



A'i{x-xi,X2,C) = Al£a-x + ^^-^{- -] (a = 3,8) (5.15) 

Att \|a; — a^il \x — X2\ J 

As in the non-relativistic QED case (3.8) the field depends on the positions of the 
quarks, and now also on their color. Since (5.15) is a solution of the EOM the action is 
stationary under local variations of A^, for any constants and unit vectors £a- However, 
variations of these parameters is a global variation which can affect the action. In fact, 

= E f^^' / + T^'^ia ■ {X, -x,) + [g^) 

a=3,8 L 

The parameter "^^^^ g A^ is multiplied by the (infinite) volume of space. This term does 
not affect bound state evolution provided it is the same for all Fock components. Hence 

A' = E (5.17) 

a=3,8 

should be a universal constant, independent of xi,X2 and the quark color C. The O (g) 
interference term is finite and was evaluated as in (2.4). It is stationary wrt. variations of 
the unit vectors £a provided £a \\ xi — X2. Choosing £a = T^*"(a;i — X2)/\T^''" {xi — X2)\ 
gives (as seen below) an attractive linear potential oc gA'^\T^''~' {xi — X2)\ between quarks 
of color C. 

The (instantaneous) action (5.16) should be stationary also wrt. variations in the ratio 
As/Ag which leaves A in (5.17) invariant. Using a lagrange multiplier A, the extremum of 
the O (g) term in (5.16) for quark color C = 1, 

Sgf = I (aI + ^A^) - 0.21 + A(A4 - A| - A|) (5.18) 

for variations of A3, Ag and A gives Ag/A^ = -v/3 and thus 

oA^ 

Sgf = ^\x,-X2\ (5.19) 

The calculation and result is the same for C = 2, whereas for C = 3 the extremum of 

Sgf = - + A(A4 -At- At) (5.20) 

is obtained for A3 = 0, Ag = A, giving Sg^^ = Sg^^ = Smt- The fact that the stationary 
value of the interference term is independent of quark color is a consequence of the color 
singlet nature of the action (5.16) and the color covariance of the EOM. 
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Having determined the parameters Aq and ta in the potential (5.15) for each Fock 
state we may now proceed to impose a stationary time dependence on the bound state, as 
already indicated in (5.8). Analogously to the non-relativistic case (3.14) we have 



^ = R(0|i ^-j^ '(/'na(0,a:;i)|£^,0) + zV'rf^' (0, 2:2) 1^,0) 

+Rm'^^p{Q,X2)i^^Mxi)[Hm - S^nt]\E,Q) = EcP^^{0;xi,X2) (5.21) 



where the energy in the field contributes —Smt as shown in (3.13). The interaction hamil- 
tonian 

Hlit) =9^2 I d'^i^fit^ ^)^a(^) T^^'i'fit, ^) (5.22) 

f=u,d 

is diagonal in color for the field (5.15) and thus consistent with the color structure (5.13) 
of the wave function. Its matrix element in the bound state equation (5.21) with C = D 
contributes (no sum on C, and neglecting terms of O (5^)), 

R(0|vg(0,a;2)VS.(0,a:i)FK0)|i?,0) = 5 J^^f^ [AI{x^) - AUx^)] xi^i^x^) 

a 

ah? 

= gY,^l\Ta^ {Xi - X2)\X{X1,X2) = ■^\xi-X2\=^Sint (5.23) 
a V 3 

Thus the interaction energy is independent of the quark color C, similarly to the instanta- 
neous action Sint- 

Using (4.6) the bound state equation (5.21) for the color singlet ud wave function 
becomes 

l^{-iVi--f + mu)x{xi,X2)-x{xi,X2)-i^{iV2-l + md) = [E-V{xi,X2)]x{xi,X2) (5.24) 

which has the same form as (3.17) for non-relativistic QED atoms. Due to the use of the 
retarded vacuum (5.2) this equation may be applied also to relativistic bound states at 
lowest order in K and to O (g) in the gauge coupling with the linear potential 

V{xi,X2) = \xi - X2\ (5.25) 

where A is a free parameter with dimension of mass. Separating the CM momentum k 
according to 

X{xi,x2) = e^^-(-i+-^)/2 xkixi - X2) (5.26) 

the bound state equation reduces to the form (2.19) given in Section 2.3. 

The bound state equation (5.24) is a rather natural generalization of the Dirac equation 
and as such has been studied before [5, 6, 7, 9]. As mentioned in Section 2.3 it has several 
intriguing properties, in particular a correct dependence of the bound state energy E on the 
center-of-mass momentum fc, and rapid oscillations of the wave function at large distances 
r between the quarks, where V{r) ^ E. 
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6. Baryon bound state equation in QCD 

The derivation of the baryon uds equation foUows the same principles as that of mesons, 
giving a specific three quark potential. I assume distinctly flavored quarks for simplicity. 
The baryon state at t = is expressed as 

3 

\E,t = 0)= I Hd^^i^^lit = 0,2/i)^fJ^(t = 0,y,)i;%{t = 0,y3)x^^BT'(?/i>y2,?/3)|0>K 
i=i 

(6.1) 

where now^ 

|0)/? = H dt\p, A) (p, A) dj\p, A)|0) (6.2) 

p,A,A 

The baryon state (6.1) is invariant under time independent gauge transformations il^'^{t, x) — >■ 
U^^' {x)il)'^' {t,x) provided the wave function is transformed as 

XABC{XI,X2,X^) ^ U"^^' {xi)U^'^' {X2)U^^' {X^)XA'B'C'{X1,X2,X^) (6.3) 

I assume that there is a gauge where the wave function has the standard color dependence 
X^^B"a"^(^i,^2,^3) = eABcx"^"^"n^i,^2,a;3) (6.4) 

As for mesons, this allows O (^g^^ instantaneous gauge fields only for a = 3, 8. The color 
ABC = 123 Fock state matrix elements of the QCD equations of motion 



f=u,d,s 

give for a = 3, 8 and u = Q 



\E,t) = 
(6.5) 



3 

-V2A0(a;) + ^0 = gY^T^^^x - x^) {a = 3,8) (6.6) 

i=i 

The commutator term oc fabc does not contribute at O (g) when the O (^g^^ fields appear 
only in the commuting elements (6, c = 3, 8) of SU(3). The solution including homogeneous 
linear terms are then to O (g) 

Alix- {x,}, ABC = 123) = Alh-x + fl ( " V^^) (6-7) 

47r 2 \ |a; — a^il \x — X2\J 

Al{x;{xi}, ABC = 123) = Klh-x + ^^ f , ^ , + ^ 



2^/?, \\x — Xi\ \x — X2\ \x — X2,\ 

This gives the instantaneous action corresponding to (5.16), 

-i E / d'^K^FH^ = E [ d'x + + O (g^) 

a a=3,8 ^ 



"This definition of the retarded vacuum is different from (5.2) for mesons. This incompatibihty needs 
to be addressed in order to treat meson-baryon interactions. 
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Similarly to the meson case the O [g^) term proportional to the (infinite) volume of space 
must be universal, implying the constraint (5.17). The O (g) interference term for the color 
component ABC = 123 is 



6 '' + 6^3 



^3 _X2\ + ^ \xi + X2- 2x-i\ + A(A4 - A^ - A^) (6.9) 



where in the second line I used stationarity of 5*/^^"^ to fix the directions of the unit vectors, 



3 II xi — X2 and £g \\ xi + X2 — 2x'^, and added the lagrange multiplier for the constraint 

'123 
int 



(5.17). The extremum of Sj^f is obtained with 



^ = V3 , , (6.10) 



giving 



Slaf = 01 ^Jxl+xl + xl-Xl■X2-X2■X3-X3■Xl (6.11) 

This expression is fully symmetric under permutations of xi ^ X2 X3, ensuring that the 

same result will be obtained for all color components of the wave function (6.4): Sj^^ = 

q213 _ = Q 
'-'int — ■ ■ ■ — >->mt ■ 

The stationarity condition for the baryon state (6.1), 

i^\E,t) = E\E,t) (6.12) 

imposes an O {g) condition on the wave function x^BC'"^(a^i5 X2,X3) which is analogous to 
(5.21) for mesons. In the interaction Hamiltonian (5.22) only the linear, O (5'^) terms in 
the gauge fields (6.7) need be considered, with the parameters £3, £s and Aa/Ag determined 
as above by the extremum of the action for each Fock state. The stationarity condition is 
diagonal in color and for the ABC = 123 color component reads 

3 3 
Yl [7°(-^V,- • 7 + mj)] x(.x,,X2, X3) + gY,Yl Ta'K{^j)x = {E + Sint)x (6-13) 

j=l j=l a=3,8 

where the interaction term on the Ihs. is 
oA^ - oA^ 

£3 • {Xi - ^2) + 1^ 4 • {Xl + X2- 2X3) = ?,S,nt{xi,X2, X3) (6.14) 

The fact that Smt given by (6.11) is a symmetric function of the quark positions Xj implies 
that the potential is the same for all color components and thus compatible with the color 
structure (6.4) of the wave function. The bound state equation for the uds baryon wave 
function x{xi-,X2,X3) is then 

3 

[7°(-^V, • 7, + m,)] x = {E- V)x (6.15) 
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where Vj = d/dxj, the Dirac matrices 7^- multiply the jth Dirac index on x' 
"^u, 1^2 = 'nT'd, TTT'S = "Is and the potential may be expressed as 



, mi = 



V{xi,X2,X3) 



3\/3 



\/ {xi - X2f + (032 - 2:3)2 + (aja 




(6.16) 



In the limit where two quarks are in the same position, e.g., X2 = 3^3, this potential coincides 
with the meson potential (5.25). It is perhaps significant that a potential for more than 
three quarks cannot be constructed analogously, as there are then more independent quark 
separations than diagonal SU(3) generators. 

7. Discussion 

I presented an approach to gauge theory bound states that uses hamiltonian time evolution 
in Coulomb gauge at lowest order in ^, i.e., in the absence of loop corrections. The ^ — )• 
limit selects gauge field configurations for which the action is stationary. At each instant 
of time, and for each Fock component of the wave function, J^(x) is determined by the 
positions of the charged constituents. The absence of loops (Born level) allows the use of 
retarded boundary conditions. This avoids the explicit appearance of higher Fock states 
arising from time ordering (Z-graphs), without affecting the energies of the bound states. 

The field equations determine J^{x) only up to a homogeneous solution, which allows 
adding a linear term to the Coulomb potential. Such a term would not be acceptable in 
the form of a fixed external potential, since it would break rotational invariance and imply 
a field extending to spatial infinity. However, J^(x) is a non-propagating field which is 
determined by the positions of the charged constituents at each instant of time. Stationarity 
of the action requires a linear potential term to be oriented along the charge separation in 
each Fock component, ensuring rotational invariance. For neutral states the coherent sum 
of the contributions to J^{x — t- 00) from the various Fock states cancels, thus avoiding a 
long-range field external to the bound state. 

A linear potential implies confinement and is thus of obvious interest for QCD. The 
present semi-classical approach does not predict its magnitude nor why such a term should 
be present in QCD but not in QED. As seen from (2.3) a linear term gives a contribution 
to the instantaneous action which is proportional the volume of space, and thus appears 
as a boundary condition related to vacuum structure which is motivated by data. Its 
derivation from theory would require methods beyond the semi-classical approximation, 
such as lattice gauge calculations or the use of Dyson-Schwinger equations [12, 13]. 

In a purely Coulombic bound state the constituent velocity scales as f ~ a. Hence 
relativistic motion u ~ 1 requires a large coupling a ^ 1, for which the use of perturbative 
methods cannot be justified. A linear potential of strength gb^ dominates, in a perturbative 
sense, gluon exchange interactions of O (5^)- This enables the study of relativistic bound 
states even for small values of the gauge coupling g. The phenomenological success of the 
Quark Model, which uses a combination of a linear and an O (a^) perturbative potential, 
indicates that the strong coupling may freeze at a sufficiently small value to justify the use 
of a perturbative expansion even for hadrons. 
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Symmetries of the lagrangian are maintained at each order of h and g. The frame 
dependence of wave functions defined at equal time of the constituents is non-trivial, since 
boosts are dynamical operators which do not commute with the hamiltonian. In the case of 
a purely linear potential the energy eigenvalues E of the bound state equation (2.19) were 
found [9] to have the correct dependence on the CM momentum k. The wave functions are 
related by a generalized Lorentz contraction (2.20), with a length scale proportional to the 
inverse canonical energy 1 / (E — V) instead of 1/E. To my knowledge this is the only case 
where an explicit relation between equal-time wave functions in different frames has been 
established, and it deserves further study. It may allow to relate spherically symmetric wave 
functions in the rest frame with the infinite momentum frame wave functions corresponding 
to quantization on the light-cone. The latter wave functions are invariant under boosts but 
defined wrt. a fixed direction in space [14]. 

The wave functions have a rapidly oscillating phase ~ exp{irV (r) / 4) at distances 
r between the constituents where the linear potential V{r) ^ E. Hence the relative 
momentum between the bound state constituents can be large, allowing support for quark 
distributions at low x in DIS, as well as Regge type quark exchange between hadrons 
scattering at high energies. This is interesting also in view of the linear Regge trajectories 
found in the spectrum [6]. 

The meson and baryon solutions found here are essentially equivalent to those in [8], 
which relied on some ad hoc assumptions. The present framework is based on a power 
expansion in h and g and thus close to a systematic derivation of bound states from the 
lagrangian. This should allow to address further issues, checking the viability of the method 
presented here. Many open questions obviously remain. A central outstanding issue is 
whether spontaneous chiral symmetry breaking can be described at the semi-classical level, 
e.g., using a chirally non-invariant vacuum state. It is also essential to demonstrate that 
loop corrections can be systematically evaluated. 
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Appendix 



A. Time development using the functional integral 



In this paper I derived bound state equations using a hamiltonian method with boundary 
conditions specified by the retarded vacuum (2.15). The time development of fermion 
states may equivalently be considered in the functional integral formulation. The boundary 
condition is then imposed through vacuum wave functionals V{ti), y{tf) at the initial and 
final times [11], 

{f\0\i)= r V[^,^]V\tf)OV{ti)exp[iS/h] (A.l) 



The Feynman propagator Sf{x — y) is obtained for O = '0(x)'0(2/) using a wave functional 
Vp with a Dirac sea of filled negative energy states i.e., the usual perturbative vacuum, 



Vpit) = exp 



exp 



d^p 1 
(27r)3 2E~p 

1 r d^p 



'^{t,p){p ■ 1 + m)il){t,p) 



(A.2) 



(2vr) 



In the second expression the wave functional was diagonalized using the positive and neg- 
ative projections 



r{t,p) = Y.Ulx^P±At,p) ; ui, = ^ 



v?{-p) 



(A.3) 



where the ± refers to positive/negative energy and A = it^ to the spin states. The minus 
sign in front of ^{t,p)'ip-^x{t,p) in (A.2) signifies that the negative energy states are 
filled. This distinction between positive and negative energy field components implies that 
Vp is non-local in coordinate space. 

The retarded propagator Sji{x — y) (2.12) is obtained using the wave functional Vr 
that has unfilled negative energy states. 



VR{t) = exp 



d-^p 
(2^ 



Yl [V'i,A(*'P)V'+,A(t,p) + ^L,A(i'P)V'-,A(i,P) 



exp 



d^xil){t,x)-i^^{t,x) 



(A.4) 



and hence is local in x. 

The functional integral method is well suited for studying the ^ — t- limit since the 
exp[iS'/^] factor in (A.l) restricts the action S to a. stationary value, thus fixing the gauge 
field at its classical value for a given charge distribution. 
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